AFCRL-TR-76-0057 

ENVIRONMENTAL  RESEARCH  PAPERS.  NO.  550 
HAES  REPORT  NO.  38 


r an 

Wm 


rNj  Rocket  Measurement  of  OH  Emission 
0^  Profiles  in  the  1.56  and  L99yu.m  Bands 


W.  F.  GRIEDER 
K.  D.  BAKER 
A.  T.  STAfR,  Jr. 
J.  C.  ULWICK 


28  January  1976 


Approved  for  public  release;  distribution  unlimited. 


OP15CAL  PHYSICS  LABORATORY  PROJECT  7670 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 

HANSCOM  AFB,  MASSACHUSETTS  01731 


AIR  FORCE  SYSTEMS  COMMAND,  USAF 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
Documentation  Center.  All  others  should  apply  to  the  National 
Technical  Information  Service. 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whm  0«a  Enl,  • O 


REPORT  DOCUMEHTATIOH  PAGE 


in;  OT.iCTTnia 


AFCRL-TR -76-0057.  A FCK 


HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


CIPIENT'S  CATALOG  HUMBER 


S.  TYRE  OF  REPORT  * PERIOD  COVERED 


S.  PERFORMING  ORO.  REPORT  NUMBER 

ERP  No.  550  ^ 


>■  CONTRACT  OR  GRANT  NUMBEA(c) 


»•  RCRPORMING  ORGANIZATION  NAME  AND  ADDRESS 


Air  Force  Cambridge  Research  Laboratories(OP 
Hanscom  AFB  STZ 

Massachusetts  01731 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ^ ] It.  REPORT  OATS  c — — 

Air  Force  Cambridge  Research  Laboratories(OpJ-^  28  Janm*p**76 ) 
Hanscom  AFB  (f  ljk.  number  of  pages  ' 

Massachusetts  01731  VjsT  73 

U.  MONITORlHO  AGENCY  NAME  A ADDRESVII  dlllttnl  from  Conlnlllnt  officl)  15.  SECURITY  CLASS,  fol  Mila  ..port; 


Unclassified 


I5a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


I IS.  DISTRIBUTION  STATEMENT  for  Oil.  RafUNtj 


Approved  for  public  release;  distribution  unlimited. 


10.  II  Igffaranu 


IS.  SUPPLEMENTARY  NOTES 

♦Utah  State  University 
Logan,  Utah  84322 


19>  KEY  WOROS  (Confirm*  on  rovoroo  side  It  n*c***ary  mut  Identity  by  block  ntmbor) 

OH  emission  r 

SWIR  backgrounds  [ s ; . 

D region  aeronomy  f L 

Infrared  " . 


r q Jip*  V 


20.  ABSTRACT  (Confirm*  on  r*v*r**  oldo  It  nocottory  md  Identity  by  block  numbor) 

" ^"Results  of  radiometer  measurements  of  infrared  Hydroxyl  (OH)  emissions 
from  five  rocket  flights  are  presented.  Two  of  the  rocket  payloads  were 
flown  from  Poker  Flat  Rocket  Range,  Alaska  and  three  were  flown  from  White 
Sands  Missile  Range,  New  Mexico.  Four  of  the  flights  were  flow*)  during 
nighttime  conditions  and  one  (at  WSMR)  was  flown  during  twilight  (x  = 99°)) 
Each  rocket  payload  wa;  instrumented  with  liquid  nitrogen  cooled'anaF"""' 
channel  radiometers  to  measure  infrared  OH  emissions  from  the  vibration- 
rotation  bands  2-0  through  5-3  (1.  5 6 (^im)  channel)  end  the  7-5  through  9-7 

0°  , SSPti  1473  «wth»n  OF . NOV  .s o.solete/  Unclassified 

y * 1 CURITY  CLASSIFICATION  OF  THU  PAGE  (»)\tn  D«l«  Ent trod) 


('  ' v~  ^ 


V /*  \ 

s>  v.' 


Unclassified 


SECURITY  CLASSIFICATION  OF  This  PAGEfWiKl  Dal*  Bnttrtd) 


/ 


20.  Abstract  (Continued) 

(1.  99(/imv)channel).  The  rocket  data  shows  that  the  principle  OH  emissions  in 
the  bands  cited  occurs  in  a thin  layer  8 to  10  km  thick  with  pe^x  emission  near 
an  altitude  of  89  km.  The  ratio  of  nighttime  1.  56\Mrryto  1.  99uim) radiance  is 
constant,  at  a value  of  1. 4,  up  to  an  altitude  of  80  km  regardless  of  latitude. 

At  twilight,  this  ratio  is  larger  at  low  altitudes  (60  km)  and  diminishes  with 
increasing  altitude.  -Estimated  nighttime  peak  volume  emission  rates  at  WSMR 
9je.4-megaphotojs  ^crn  ~.se£'i..ttm  "*,in  the  1.  56t^m  band  and  0.7  megaphotons 
(cm-3  sec"l  Mm“lHn  the  1.  99;pm;b4'nd.  _ \ J 

' ’ i 

i 


A 


) 


i 


j 

I 

| 

I 


Unclassified 

SECURITY  CLASS!  CATION  OF  THIS  FAOEfWi*"  Out  tnfrU) 


Foreword 

A High  Altitude  Effects  Simulation  (HAES)  program  is  being  conducted  by  the 
Defense  Nuclear  Agency.  Under  HAES,  several  related  but  distinct  measurement 
programs  have  been  or  will  be  conducted  to  provide  information  for  development 
and  test  of  predictive  computer  codes  which  are  used  to  assess  and  evaluate  the 
operation  of  critical  DOD  radar  and  optical  infrared  systems  in  nuclear  disturbed 
environments.  The  Air  Force  Cambridge  Research  Laboratories  work  under  the 
HAES  program  will  be  reported  under  the  HAES  Series,  AFCRL  Environmental 
Research  Papers.  The  reports  will  also  identify  the  measurement  program  (for 
example,  ICECAP,  EXCEDE,  SPIRE,  etc. ) so  that  related  reports  and  results  can 
be  correlated  and  utilized.  One  of  the  measurement  programs,  ICECAP  (Infrared 
Chemistry  Experiments  — Coordinated  Auroral  Program),  is  conducted  in  the  polar 
disturbed  atmosphere.  The  specific  objective  of  ICECAP  is  to  execute  coordinated 
field  measurement  programs  to  investigate  ionization  and  excitation  mechanisms 
and  chemical  processes  leading  to  both  short  and  long  wavelength  infrared  (SWIR, 

1 to  6 /urn  and  LWIR,  6 to  30  nm)  emissions  in  an  auroral  display.  The  measure- 
ments are  coordinated  utilizing  rocket,  balloon,  aircraft  and  ground-based  obser- 
vational platforms  where  appropriate.  The  program  has  at  present  three  distinct 
phases:  (1)  Observation  of  infrared  and  related  emissions  along  with  production 
sources  and  ionospheric  measurements  in  bright  auroral  arcs;  (2)  Observations  of 
infrared  aurorally  associated  emissions  during  daylight;  and  (3)  Observations  of 
hydroxyl  (OH)  auroral  enhancements. 

The  report  presents  rocket  measurements  of  infrared  emissions  in  the  vibra- 
tion-rotation bands  of  the  hydroxyl  (OH)  molecule  in  both  mid-latitude  and  auroral 
zones.  Particular  emphasis  is  given  to  data  reduction  procedures,  correction  of 
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data  for  rocket  aspect,  and  computation  of  OH  volume  emission  rates.  A con- 
certed effort  is  made  to  present  all  pertinent  information  necessary  for  other 
analysts  to  utilize  the  measured  data  for  individual  studies. 

A.  T.  STAIR,  JR.  JAMES  C.  ULWICK 
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Rocket  Measurement  of  OH  Emission 
Profiles  in  the  1.56  and  1.99yu.m  Bands 


1.  INTRODUCTION 

Under  the  sponsorship  of  the  Defense  Nuclear  Agency  several  rocket  measure- 
ment programs  (ICECAP  Programs)  have  been  accomplished  in  the  past  few  years 
to  investigate  infrared  emissions  in  the  upper  atmosphere,  particularly  under  dis- 
turbed conditions  such  as  during  aurora.  Part  of  this  effort  has  been  dedicated  to 
studies  of  emissions  from  the  hydroxyl  (OH)  molecule  using  liquid  nitrogen  cooled 
radiometers  onboard  small  rockets.  Rocket  measurements  are  employed  in  order 
to  determine  the  altitude  distributions  of  infrared  emissions  from  which  the 
mechanisms  controlling  the  degree  of  OH  excitation  in  the  atmosphere  are  deduced. 
Results  from  the  initial  rocket  experiment  under  the  ICECAP  program  have  been 
reported  previously  (Rogers  et  al1).  This  report  presents  the  results  from  five 
additional  rocket  measurements  of  OH  emission  profiles  using  similar  payloads. 
Two  of  these  payloads  were  flown  from  Poker  Flat  Rocket  Range  (PFRR),  Alaska 
during  the  ICECAP  73  program  (March,  April  1973)  and  three  were  flown  from 
White  Sands  Missile  Range  (WSMR),  New  Mexico  in  October  1973  during  a joint 
Army/Air  Force  program  called  "Midlatitude  Twilight  D-Region  Studies. " 


(Received  for  publication  23  January  1976) 

1,  Rogers,  J.  W. , Murphy,  R.E. , Stair,  Jr.,  A.T. , and  Ulwick,  J. C.  (1973) 
Rocket-borne  radiometric  measurements  of  OH  in  the  auroral  zone, 

J.  Geophys.  Res.  78:7023. 
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The  technique  employed  in  these  experiments  was  to  measure  the  infrared 
emissions  in  the  1.  56  and  1.  99-^m  bands  with  a dual  channel  radiometer  payload 
as  the  rocket  traversed  the  CH  emitting  layer  (70  to  95  km).  The  bandpasses  of 
the  two  channels  were  selected  to  encompass  the  2-0  through  5-3,  and  the  7-5 
through  9-7  OH  vibration-rotation  bands,  respectively. 

These  particular  optical  bands  were  chosen  in  order  to  give  information  about 
the  excitation  processes  responsible  for  the  OH  emissions.  In  particular  the  ozone 
process 

Og  + H — • OH*  + 02 

can  energetically  excite  up  to  and  including  the  ninth  vibrational  level  whereas 
other  possible  processes  cannot  excite  as  high  levels.  For  example,  the  mechan- 
ism 

O + H02  - OH*  + Oz 

can  only  excite  v * 6 and  lower  levels.  As  a result,  measurements  of  the  emission 
in  the  two  bands  can  give  a relative  measure  of  the  completing  mechanisms.  This 
report  documents  the  rocket  data  for  these  experiments  and  the  data  processing 
techniques  used  to  determine  the  OH  volume  emission  rates.  In  addition,  some 
supporting  measurements  from  ground-based  instruments  conducted  during  the 
rocket  launches  are  presented. 


2.  INSTRUMENTATION 

2.1  Rocket  Payload  Configuration 

Determination  of  the  altitude  distribution  of  the  OH  emission  can  be  accom- 
plished by  measurement  of  the  integrated  overhead  infrared  emissions  through  the 
altitude  region  from  70  to  95  km.  The  rocket/payload  combination  used  for  this 
purpose  was  a single-stage  solid  fuelAstrobee  D equipped  with  cryogenically  cooled 
dual-channel  radiometers. 

The  rocket  is  a 6-in.  diameter,  fin-stabilized  vehicle  which  provides  adequate 
acceleration  to  offset  surface  wind  effects,  yet  exhibits  a relatively  low  sustained 
acceleration  to  minimize  environmental  problems.  It  is  capable  of  carrying  a pay- 
load  weight  of  30  lb  to  approximately  95  km  in  140  sec  (sea  level  conditions). 
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The  primary  sensor  used  for  the  OH  measurements  was  a Utah  State  Uni- 

2 

versity  dual-channel  radiometer  (Jensen  et  al  ).  A schematic  of  the  radiometer 
is  shown  in  Figure  1.  It  consists  of  an  optical  subsection  containing  Indium 
Antimonide  (InSb)  detectors,  collecting  optics,  and  interference  filters  in  a cryo- 
genic dewar  cooied  to  near  liquid  nitrogen  temperature  (77°K).  The  components 
provide  tv/o  independent  optical  channels,  which  utilize  a common  optical  chopper 
to  modulate  the  incident  radiation.  The  system  has  an  ejectable  cold  cover  to  keep 
the  optical  system  cold  and  yet  protect  it  from  frosting  until  a suitable  altitude 
(approximately  50  km)  where  the  cover  is  ejected  along  with  the  payload  nose  tip, 
thereby  exposing  the  radiometers  (see  Figure  2).  The  optical  bandpass  of  each 
channel  is  determined  by  the  interference  filters,  nominally  Aq  = i.  5642, 

AA  = 0.  2057  pm  for  the  low  OH  vibrational  levels  and  A^  = 1.  9850,  AA  = 0.  2027  pm 


Figure  1.  Schematic  of  Liquid  Nitrogen  Cooled  OH  Radiometer  (Jensen  et  al2) 


2.  Jensen,  L.  L, , Kemp,  J.  C. , and  Bell,  R.J.  (1972)  Small  Rocket  Instrumenta- 
tion  for  Measurement  of  Infrared  Emissions.  Astrobee  D 30.  205-3  and 
Astrobee  D 30.  205-4,  Sci.  Rept. , No.  3,  AFCRL  72-0691,  Contract  No. 
F19628-70-C-0302,  Utah  State  University,  Logan. 
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Ejected 

Nosetip 


TENTATIVE  PAYLOAD 

ASTROBEE  D'«  A30.205-3,  A30.203-4 
ICECAP  73 


Figure  2.  Astrobee  D Rocket  Payload  Configuration 
Showing  Nose  Tip  Ejection  Technique  and  Nominal  Fields 
of  V iew 


for  the  high  OH  vibrational  levels.  The  nominal  system  spectral  response  for  each 
channel  is  shown  in  Figures  3 and  4 (Wyatt  and  Kemp  ).  A tabulation  of  the  sys- 
tem spectral  response  data  is  included  in  Appendix  A.  The  system  wavelength  pass- 
bands  were  computed  by  integrating  the  system  spectral  response  and  normalizing 
the  result  to  an  ideal  filter  of  square  response.  For  convenience,  some  references 
in  this  report  to  these  passbands  and  center  wavelengths  will  be  shortened  to 
AA  = 0.21  pm  for  the  A.  = 1.56-pm  channel  and  AA  = 0.20  pm  for  the  1.  99-pm 
channel. 


3. 


Wyatt,  C.  L. , and  Kemp,  J.  C.  ( 
NR-3B-8,  NR-3B-9,  NR-3B-10, 


) Calibration  of  SWIR  Radiometers,  Model 
NR-3B-11,  NR-3B-12. 
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RELATIVE  RESPONSE 


Wavelength  Channel.  Center  wave- 
length X = 1.  985  pm,  equivalent 
spectrafpassband  AX  = o.  2027  pm) 


Figure  4.  Radiometer  Relative  System 
Spectral  Response  for  Short  Wavelength 
Channel.  Center  wavelength  X = 

1,  5642  pm,  equivalent  spectral  passband 
AX  o.  2057  pm) 


The  radiometer  optics  consists  of  a simple  f/0.  58,  high-transmission,  coated, 
silicon  lens  and  establishes  the  field-of-view  of  approximately  5°  full  angle  (see 
Appendix  A for  measured  fields  of  view).  A light-emitting  diode  (LED)  within  the 
optical  section  is  periodically  activated  to  stimulate  the  detector  and  provide  an 
inflight  indication  that  the  system  is  responding  normally.  Because  of  the  variable 
nature  of  LED  devices,  the  inflight  radiometer  signals  resulting  from  them  cannot 
be  used  for  absolute  calibration  of  the  radiometer  but  serve  only  as  a genera]  indi- 
cation of  system  performance.  Special  techniques  were  required  to  calibrate  the 
radiometers  and  have  been  reported  previously  (Wyatt  and  Kemp  );  however,  the 
pertinent  results  of  these  calibrations  are  included  in  this  report. 

The  electrical  bandwidths  of  the  radiometer  channels  are  1 Hz  and  consequently 
the  instrument  integration  time  is  of  the  order  of  a second. 

In  addition  to  the  radiometer,  each  payload  was  equipped  with  an  aspect  system 
consisting  of  a magnetometer  (Schonstedt  Instrument  Company)  and  an  infrared 
horizon  sensor.  For  the  auroral  zone  payloads,  a photometer  measuring  N2  emis- 
sion at  3914A  was  also  included  to  assess  auroral  conditions  prevailing  during  the 
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OH  measurements.  Temperature  sensors  were  included  in  all  payloads  to  moni- 
tor the  environment  in  critical  areas  of  the  payload  and  to  assess  instrument  per- 
formance. 

Standard  frequency  deviation  multiplexing  telemetry  (FM/FM)  techniques  were 
used  to  transmit  the  sensor  data  to  the  ground.  The  analog  voltage  signals  (0  to 
+5V)  from  the  rocket  sensors  were  used  to  frequency  modulate  standard  IRIG  sub- 
carrier oscillators  (IRIG  channels  12  through  16)  - each  data  output  was  assigned  a 
separate  subcarrier  frequency.  The  modulated  subcarriers  were  then  multiplexed 
into  a single  composite  signal  which  was  used  to  frequency  modulate  an  S-band 
telemetry  transmitter.  The  ground  telemetry  stations  received  the  S-band  carrier 
and  detected  the  FM  composite  signal  containing  the  measured  information.  The 
composite  signal  was  direct  recorded  on  magnetic  tape  at  60  in.  per  second  along 
with  IRIG  B format  Universal  Time  Code.  The  recorded  single  composite  signal 
containing  the  data  from  all  the  sensors  on  the  rocket  is  the  starting  point  for  proc- 
essing to  retrieve  the  measured  information. 

2.2  Ground-Based  Instrumentation 

A variety  of  ground-based  optical/infrared  and  ionospheric  sounding  instrumen- 
tation was  operated  in  support  of  the  rocket  programs  both  at  PFRR  and  WSMR,  At 
both  locations,  the  AFCRL  Mobilab  and  the  Utah  State  University  Mobile  Observatory 
(ARGUS)  were  operated.  Instrumentation  in  these  facilities  included  interferometer- 
spectrometers,  photometers  (fixed  and  spatially  scanning),  all  sky  cameras, 

30-MHz  riometer  and  several  near-infrared  radiometers.  At  PFRR,  the  University 
of  Alaska  provided  additional  optical  support  as  well  as  magnetometer  coverage. 

A complete  description  of  these  systems  and  interpretation  of  results  obtained  is 
beyond  the  scope  of  this  report;  however,  some  data  are  included  in  the  text  where 
relevant  to  define  the  geophysical  conditions  prevailing  during  the  rocket  flights. 


3.  ROCKET  FLIGHT  SUMMARY 

Table  1 summarizes  the  rocket  flights  flown  during  the  OH  measurement  pro- 
grams at  PFRR  and  WSMR.  Other  rocket  experiments  were  also  conducted  during 
these  field  programs,  but  are  not  reported  here. 

3.1  Rocket  Flight  Characteristics 

An  understanding  of  rocket  flight  characteristics,  particularly  the  instrument 
viewing  aspect,  is  critical  to  correct  interpretation  of  the  infrared  measurements. 
Of  particular  importance  is  the  Increase  in  signal  commonly  referred  to  as 
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Table  1.  Summary  of  OH  Rocket  Launches 


Launch  Site^ 

Rocket 

Number 

Launch  Data/ Tim  e^ 
(UT) 

Launch  Conditions 

PFRR 

A30. 205-5 

080.  1011:00.  000 
(21  Mar  1973) 

Nighttime 

Normal  OH  Backgroun  , 

PFRR 

A30.  205-6 

096.  0845:00.000 
(6  April  1973) 

Nighttime 

Quiet  OH  Background 

WSMR 

AO30.  311-1 

276.0040:00.  949 
(4  Oct  1973) 

Sunset 

WSMR 

AO30.  311-2 

276.  0127:00.818 
(4  Oct  1973) 

Twilight  OH 

WSMR 

AO30.  311-3 

276. 0500:00.716 
(4  Oct  1973) 

Nighttime  Quiet  OH 

Rocket 

Number 

Solar  Zenith 
Angle 
(deg) 

Apogee 

(km) 

n 

Altitudes  of 
Usable  Data 
(km) 

A30.  205-5 

114.7 

78 

55-78  Ascent 
78-25  Descent 

A30.  205-6 

107.5 

78 

53-78  Ascent 

AO30.  311-1 

56 

52«> 

None<3) 

AO30.  311-2 

99 

102 

52-78  Ascent 

AO30.  311-3 

141 

106 

55-106  Ascent 
106-74  Descent 

(1)  PFRR  - Poker  Flat  Rocket  Range,  Alaska;  WSMR  — White  Sands  Missile  Range, 

New  Mexico. 

(2)  Date/Time  - Julian  day  hours,  minutes:  seconds,  milliseconds  UT. 

(3)  Early  tip  ejection  resulted  in  loss  of  data. 


van  Rhijn  effect  (Chamberlain  ),  when  an  emitting  region  is  viewed  obliquely  rather 
than  vertically.  The  optical  axis  of  the  radiometer  was  aligned  with  the  rocket  axis 
in  each  case  so  that  any  inclination  of  the  rocket  with  respect  to  the  vertical,  allows 
the  radiometer  to  view  a greater  optical  path  through  an  emitting  layer  resulting  in 
a larger  measured  signal, 


4.  Chamberlain,  J.  W.  (1961)  Physics  of  the  Aurora  and  Airglow,  Academic  Press, 
New  York  and  London. 
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A scenario  of  a typical  Astrobee  D rocket  flight  is  illustrated  in  Figure  5. 

The  rocket  is  launched  at  an  elevation  angle  of  about  85°.  At  about  20  sec  after 
liftoff,  the  motor  propellant  burns  out.  At  this  time,  the  rocket  has  attained  a 
spin  rate  of  about  8 rps  due  to  the  effect  of  the  rocket  fins.  In  addition  to  the 
rocket  spin  (roll),  the  vehicle  usually  exhibits  a coning  (precession)  motion  with 
a period  considerably  longer  than  the  spin  period,  typically  20  to  30  sec.  When 
the  rocket  has  achieved  an  altitude  of  about  55  km,  the  atmospheric  drag  is  small 
and  the  vehicle  is  essentially  in  free  ballistic  flight,  allowing  the  payload  tip  to  be 
safely  ejected  permitting  the  radiometers  to  view  the  incident  infrared  radiation. 
The  ejection  of  the  tip  usually  aggrevates  the  coning  motion.  The  vehicle  con- 
tinues this  motion  through  apogee  and  for  a time  on  descent  while  the  vehicle  is 
falling  tail  first.  Finally,  around  70  km,  the  effect  of  atmospheric  drag  perturbs 
the  vehicle  motion:  it  begins  to  tip  over,  oscillates  for  a time,  and  then  reenters 
nos*,  first.  Thus,  radiometer  airglow  measurements  are  possible  from  tip  ejec- 
tion (55  km)  to  vehicle  tip  over  (~70  km)  when  the  instruments  view  the  earth. 


Figure  5.  Schematic  of  Typical  Radiometer  Rocket  Flight  Characteristics 
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C.2  I.uuni'h  (’midi lions  nnd  Payload  Performance 

The  two  OH  rockets  flown  at  PFRR  in  March  and  April  1973  (A30.  205-5  and 
A30.  205-6)  underperformed  and  achieved  apogee  altitudes  less  than  expected.  How- 
ever, good  data  were  measured  from  tip-eject  altitude  to  apogee  and  luring  descent 
as  indicated  in  Table  1.  Supporting  ground-based  radiometric  measurements 
accomplished  during  the  night  of  2i  March  1973  indicated  a generally  normal  auroral 
zone  OH  background  of  approximately  200  kilorayleighs  in  the  wavelength  region 
A^  = l.  573  pm  and  AA  = 0.  212  pm  at  the  time  of  rocket  A30.  205-5  flight  (see  Fig- 
ure 6)  (Huppi,  5 personal  communication).  Some  auroral  activity  had  been  present 
earlier  in  the  evening.  During  the  night  of  6 April  1973  (A30.  205-6)  clouds  inter- 
ferred  with  the  ground  measurements  (see  Figure  7),  but  at  the  time  of  rocket 
launch  there  were  enough  breaks  in  the  cloud  cover  to  ascertain  that  the  OH  level 

5 

(m  the  same  spectral  band  mentioned  above)  was  reasonably  quiet  (Huppi  [personal 
communication]  estimates  a level  of  the  order  150  kilorayleighs).  During  both 
flights,  the  onboard  39 14 -A  photometers  indicate  stable  low  levels  near  the  minimum 
detectable  of  the  instruments  (~200  R). 

The  three  rockets  flown  at  WSMR  were  timed  to  study  altitude  profiles  of  OH 
emissions  associated  with  twilight-night  transition  phenomena.  Figures  8 and  9 

a 

show  the  OH  levels  measured  from  the  ground  (Huppi,  personal  communication) 
during  the  night  of  the  rocket  flights  with  radiometers  in  the  bandpasses  A^  = i.  683, 
AA  * 0.  071  and  A^  = i.  573,  &A  s 0.  212  micrometers.  Also  indicated  on  the  figures 
are  the  times  of  rocket  launches  relative  to  the  OH  levels. 

Rocket  AO30.  311-1  was  launched  prior  to  sunset  (solar  zenith  angle  was  approx- 
imately 86  degrees).  Due  to  early  tip  ejection,  this  ro:ket  underperformed  and  the 
flight  instrument  heated  up  resulting  in  a loss  of  all  mesospheric  data.  The  second 
rocket  (AO30.311-2)  was  launched  at  the  time  of  minimum  OH  emission  when  the 
OH  level  was  approximately  40  kilorayleighs  in  the  1.  683-jum  band  and  200  kilo- 
rayleighs in  the  1.  573-pm  band  (see  Figures  8 and  9).  The  solar  zenith  angle  was 
about  99°  at  this  time.  The  rocket  performed  satisfactorily  and  data  were  obtained 
from  tip  ejection  altitude  to  78  km  where  direct  sunlight  on  the  instrument  prohibited 
measurements. 

A slight  enhancement  in  the  OH  level  (ground  based  radiometer  1.  683-pm  chan- 
nel) began  at  approximately  0130  UT  and  peaked  at  a level  of  about  55  kilorayleighs 
at  0220  UT.  This  enhancement  is  typical  of  the  transition  period  and  is  also  dis- 
cernable  in  the  1.  573-pm  channel.  The  third  rocket  (AO30.  311-3)  was  launched  at 
0500  UT  when  the  OH  level  was  dropping  (44  kilorayleighs  in  the  1.  683-pm  band  and 
190  kilorayleighs  in  the  1.  573-pm  band).  This  rocket  and  payload  performed 


5,  Huppi,  R.  J.  (1973)  Private  communication. 
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Figure  6.  Ground  Based  Zenith  Radiometer  Measurement,  Poker  Flat,  Alaak; 
21  March  1973  (UT).  A = 1.573  nm.  AA  - 0,212  pm,  F.O.V.  = 2,8°  (HuppP) 
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Figure  7.  Ground  Based  Zenith  Radiometer  Measurement,  Poker,  Alaska,  6 April 
1973  (UT),  X = 1,  573  pm.  AX  = 0,212  pm,  F.O.V.  = 2.8°  (Huppi5).  Launch  of 
rocket  A30. 2lfe-6  is  indicated  at  0845  UT 


SOLAR  ZENITH  ANGLE 


Figure  8.  Ground  Based  Zenith  Radiometer  Measurement,  WSMR,  New 
Mexico,  3-4  October  1973  (UT).  X = 1.683  jum,  AX  * 0.07  pm,  F.  O.  V. 
2.8°(Huppi5)  ° 
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Figure  9.  Ground  Based  Zenith  Radiometer  Measurement,  WSMR,  New 
Mexico,  3-1  October  1973  (UT).  X = 1.  573  pm,  AX  *0.212  pm,  F.O.  V.  = 
2.  8°  (Huppi5)  0 
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nominally  and  good  data  were  obtained  from  54  km  through  apogee  (108  km)  and 
during  descent  to  74  km  where  the  rocket  tipped  over  and  reentered. 

Ground  based  riometer  (30  mHz)  data  taken  during  the  period  2 October  through 
4 October  1973  at  WSMR  indicated  that  no  detectable  ionospheric  absorption  events 
occurred  preceding  or  during  the  rocket  flights. 


4.  ROCKET  DATA  PROCESSING 
4.1  Digitization  of  Data 

The  recorded  composite  telemetry  signal  for  each  flight  was  discriminated  to 
separate  individual  sensor  channels  and  reduce  the  frequency  modulated  subcarrier 
data  to  analog  output  voltages.  Each  voltage  output  was  digitized  at  a rate  of  1660 
points/sec,  correlated  with  time  after  launch  and  coded  in  a format  compatible 
with  the  AFCRL  CDC  6600  digital  computer.  This  computer  was  used  for  all  data 
processing  and  plotting. 

Figures  10  through  17  show  the  computer  plots  of  telemetered  sensor  volts 
for  each  of  the  radiometer  channels  (denoted  by  IRIG  channel)  on  each  rocket  as  a 
function  of  time  after  launch  in  seconds.  These  represent  the  fundamental  data 
from  which  all  further  processing  was  accomplished. 

In  general,  the  data  for  each  of  the  rocket  flights  are  similar.  Prior  to  the 
nose  tip  ejection  (approximately  50  sec)  the  radiometers  view  the  cold  cover  and 
the  voltage  is  near  zero  (explanation  of  the  slightly  negative  voltages  in  some  chan- 
nels will  be  discussed  later).  After  the  tip  (and  cold  cover)  are  ejected,  the  instru- 
ments view  the  emitting  region  above  and  the  voltage  levels  show  a marked  increase. 
The  initial  large  scale  modulation  of  data  thereafter  is  due  to  coning  of  the  rocket 
which  produces  the  van  Rhijn  enhancement.  The  periodic  spikes  in  the  data  labeled 
"inflight  calibration"  result  from  activation  of  the  light  emitting  diode  for  systems 
check. 

Referring  to  Figures  16  and  17,  after  a period  of  modulated  high  voltage  levels, 
the  signal  decays  to  near  zero  because  the  rocket  (radiometer)  has  passed  through 
the  emission  region.  The  low  signal  remains  through  apogee  (around  140  sec)  and 
then  increases  again  as  the  rocket  descends  (tail  first)  through  the  emitting  regions 
once  more  (see  Figure  5).  After  220  sec,  the  rocket  has  descended  to  a low  enough 
altitude  where  the  fins  are  affected  by  air  drag  and  the  rocket  tips  over.  During 
this  period,  the  radiometers  view  first  the  horizon  and  then  the  earth.  Once  the 
rocket  tips  over  and  begins  reentry,  the  radiometer  optics  become  warm  and  satu- 
rate both  channels. 

The  absence  of  any  marki  ' signal  decay  after  tip  ejection  in  Figures  10,  11, 

12  and  13  indicates  that  these  sockets  did  not  pass  through  the  main  layer  of  OH 
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A30.205-5  IRIG  14 


Figure  10.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.5642 -jcn  Rocket  Radiometer  Channel.  Rocket  A30.  205-5  Poker  Flat, 
Alaska,  >1  March  1973  (UT) 
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Figure  11.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.  985-Mm  Rocket  Radiometer  Channel.  Rocket  A30.  205-5,  Poker  Flat, 
Alaska,  21  March  1973  (UT) 
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Figure  14.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.  5642-nm  Rocket  Radiometer  Channel.  Rocket  AO30.  311-2,  WSMR,  New 
Mexico,  3 October  1973  (UT) 


imu 


PAYLOAD  ENTERS 
•'  SUNLIGHT 

DATA 


1 i i 

WAVELENGTH 
Xo  * l.56/i  m 
AX'*  .21/1.  m 


./ 

* • •• 
•.  * JV 


i*. 


■TIP  EJECTION 


6.00 


5.00 


4.00 


<0  3.00 

H 

-I 

§ 2.00 


1.00 


0.00 
-1.00 

40  0 80.0  120  0 160  0 200  0 240.0  280.0  320.0 

TIME 

Figure  15.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.  985-^m  Rocket  Radiometer  Channel.  Rocket  AO30.  311-2,  WSMR,  New 
Mexico,  3 October  1973  (UT) 
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Figure  16.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.  5642 -(im  Rocket  Radiometer  Channel.  Rocket  A 030.  311-3,  WSMR,  New 
Mexico,  3 October  1973  (UT) 
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Figure  17.  Measured  Telemetry  Volts  Versus  Flight  Time  in  Seconds  From 
1.  986-;jm  Rocket  Radiometer  Channel.  Rocket  AO30,  311-3,  WSMR,  New 
Mexico,  3 October  1973  (UT) 
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emissions.  In  Figures  14  and  15,  the  appearance  of  signal  saturation  around 
70  sec  was  caused  by  the  incidence  of  solar  radiation  on  the  instrument  as  dis- 
cussed previously. 

The  fine  structure  in  all  the  data  is  attributed  to  instrument,  telemetry  and 
recording  system  noise  rather  than  rocket  spin  since  the  instrument  inherent  long 
integration  time  is  much  longer  than  the  rocket  spin  period. 

4.2  Measured  Altitude  Emission  Proliles 

The  voltage-time  data  in  Figures  10  through  17  were  converted  to  kilorayleigh 
radiance  profiles  using  trajectory  data  and  the  following  calibration  equation: 

Ft  = rV  + C (1) 

a 

where 

Rq  = radiance  measured  at  zenith  angle  a (kilorayleighs) 

r = inverse  responsivity  (kilorayleighs/ volt) 

V = telemetry  signal  (volts) 

C = instrument  offset  constant  (kilorayleigh) 

-2  - 1 

The  radiance  (L  ) in  watts-cm  -ster  can  be  approximated  using  the  following 
Q fi 

relationship  (Baker  ): 


X 10 


-10 


(2) 


where 

= l.  9850  fjm  for  the  high  wavelength  channel  and  A = 1,  5642  \xm  for  the 
low  wavelength  channel. 

The  appropriate  constants  used  in  Eq.  (1)  for  each  radiometer  channel  are 

O 

given  in  Table  2 (Wyatt  et  al  ).  The  offset  constant  (C)  is  a feature  of  the  differen- 
tial amplifiers  in  the  radiometer  electronics.  The  drift  of  this  offset  is  a common 
problem  with  DC  amplifiers  unless  elaborate  means  are  used  to  compensate  for  it. 
However,  as  long  as  the  product  rV  is  large  compared  to  0 in  Eq.  (1),  errors  in 
computing  radiance  caused  by  drift  in  C are  minimal.  For  this  reason  the  high 
gain  channels  (characterized  by  highest  voltages)  as  opposed  to  low  gain  channels, 
were  used  where  possible  in  reducing  the  radiometer  data.  In  all  cases,  the  offset 


6.  Baker,  D.  J.  (1974)  Rayleigh,  the  unit  of  light  radiance,  Applied  Optics  13:2160. 
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Tabic  2.  Radiometer  Calibration  Constants 


[ — — , 

X 

0 

Center 

Wavelength 

Voltage  Limits 

Rocket  Number 

IRIG 

(microns) 

(volts) 

A 30. 205-5 

16 

1.  985 

-.20  <V  <2.71 
2,71  < V < 5.  0 

14 

1.  5624 

-.45  < V < 2.  64 
2.  64  < V < 4.  80 

A30.  205-6 

16 

1*  985 

-.20  < V < 2.725 
2.725  < V < 4.  85 

14 

1.  5642 

-. 15  < V < 2.724 
2.724  < V < 5.0 

15 

1.985 

0 < V < 5.  0 

AO30.  311-2 

14 

1.  5642 

0 < V < 5.  0 

15 

1.985 

-.20  < V < 5.5 

AO30.  311-3 

16 

1.  5642 

-.35  < V < 5.5 

N 

C 

Inverse  Responsibility 

Offset  Constant 

Rocket  Number 

(kilorayleighs/volt) 

(kilorayleighs) 

56.  86 

+ 11.37 

A 30. 205-5 

627.24 

-1546.28 

- 

70.  86 

+30,47 

828.  5 

-1985.3 

57.24 

+1.  14 

A30.  205-6 

508. 78 

-1222.  1 

36.  86 

+5.  52 

413.76 

-1031.  9 

24.  19 

+7.48 

AO30.  311-2 

216.43 

-9.  93 

68.  59 

12.4 

AO30.  311-3 

34.79 

+12. 18 
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calibration  constants  were  adjusted  to  agree  with  the  voltage  level  measured  just 
prior  to  tip  ejection. 

Radiance  plots  of  the  data  as  a function  of  rocket  altitude  are  shown  in  Figures 
18,  19,  20  and  21  for  each  of  the  rockets  (in-flight  calibrations,  pre  tip-ejection 
data,  and  obvious  spurious  noise  have  been  omitted  for  clarity).  Magnetometer 
pitch  angle  data  (explained  in  the  next  paragraph)  are  also  shown  on  these  plots. 

4.3  Magnetometer  Aspect  Data  Reduction 

Each  rocket  payload  contained  a magnetic  aspect  sensor  (magnetometers) 
oriented  with  its  sensing  axis  perpendicular  to  the  rocket  axis  (see  Figure  2).  The 
rocket  magnetometer  data  alone  can  provide  only  rocket  cone  angle,  axis  of  pre- 
cession relative  to  the  earth's  magnetic  field  and  spin  and  precession  rates. 
Although  it  is  not  possible  to  completely  define  roc&et  aspect  with  magnetometer 
data  alone,  a reasonable  approximation  can  be  made  in  cases  where  "well  behaved" 

7 

flight  motion  can  be  assumed  (Marcou  ) and  additional  aspect  parameters  are 
available,  such  as  knowledge  of  the  functional  variation  of  other  measured  data  with 
aspect.  The  measured  OH  emission  altitude  profiles,  for  instance,  provide  some 
insight  into  rocket  aspect,  viz. , van  Rhijn  effect. 

The  basic  equations  relating  magnetometer  output  voltage  and  aspect  param- 
eters are  derived  in  Appendix  C.  The  pertinent  equations  are  repeated  here  in 
order  to  explain  the  data  processing  employed. 

The  instantaneous  magnetometer  output  voltage  is  related  to  both  the  angle 
between  the  earth's  magnetic  field  and  the  sensing  axis  of  the  magnetometer  and  ‘he 
rocket  spin  frequency  by 

B 

v * vave  + lf'sin  y cos  wst  (3) 

where 

v = instantaneous  voltage  output  from  the  magnetometer 

Y = angle  between  earth's  magnetic  ueld  vector  and  rocket  axis 

u = rocket  spin  radian  frequency 
s 

Bq  * total  intensity  of  the  earth’s  magnetic  field  (milligauss) 
k = magnetometer  calibration  factor  (milligauss/volt) 

7.  Marcou,  R.J.,  and  Sullivan,  B.J.  f!967)  Aspects  for  a Rocket  From  Magneto- 
meter Data,  AFCRL  67-0424,  Air  Force  Cambridge  Research  Laboratories, 
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Figure  18.  Altitude  Profile  of  Radiance  (kilorayleighs)  in  the  A 
1.  5642  u> m (TRIG  14)  and  AQ  = 1.985  Mm  (IRIG  16)  Bands,  Rocker 
A 30, 203-5.  Data  not  corrected  for  rocket  aspect 
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Figure  19.  Altitude  Profile  of  Radiance  (kilorayleighs)  ir.  the  A0  * 
1.  5642  nm  (IRIG  14)  and  AQ  = l.  085  um  (IRIG  16)  Bands,  Rocket 
A30.  205-6.  Data  rot  corrected  for  rocket  aspect 
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Figure  20.  Altitude  Profile  of  Radiance  (kilorayleighs)  in  the  XQ  * 
1.  5642  (im  (TRIG  14)  and  XQ  a i.  985  pm  (IRIG  15)  Bands,  Rocket 
AO30.  311-2.  Data  not  corrected  for  rocket  aspect 
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Figure  21.  Altitude  Profile  of  Radiance  (kilorayleighs)  in  the  X0 
1.  5642  pm  (IRIG  16)  and  XQ  ■ j.  985  pm  (IRIG  15)  Bands,  Rocket 
AO30.  311-3.  Data  not  corrected  for  rocket  aspect 


t 


= flight  time  (seconds) 

= average  magnetometer  output  voltage. 


ave 


Equation  (3)  is  an  amplitude  modulated  waveform  where  the  carrier  frequency 

is  u (due  to  rocket  spin)  and  the  amplitude  modulating  function  is  sin  7 (due  to 
s 

rocket  coning).  As  shown  in  Appendix  C,  as  long  as  7 varies  slowly  compared  to 


tuC  o^ui  w uiCH  -t  tgui  e vo/ 


sin  7=^  (vu  -vave)  (4) 

o 

where 

vy  = upper  instantaneous  envelope  voltage  from  the  magnetometer. 

The  first  step  in  reducing  the  magnetometer  data  was  to  evaluate  the  average 
envelope  voltage  (vave)  from  the  measured  data.  This  was  computed  by  simply 
averaging  the  instantaneous  maximum  (upper  envelope)  and  minimum  (lower  enve- 
lope) output  voltage  from  the  magnetometer  data  as  a function  of  flight  time.  The 
value  of  vave  was  found  to  be  reasonably  constant  for  a particular  rocket  flight. 
With  the  value  of  v&ve  established,  7 was  computed  as  a function  of  flight  time 
using  Eq.  (4).  The  value  of  the  calibration  constant  k was  supplied  by  the  manu- 
facturer (k  = 250  milligau8s/volt)  and  the  value  of  the  earth's  local  magnetic  field 
intensity  (Bq)  at  ground  level  was  used  for  rockets  flown  at  each  launch  site  (for 
WSMR,  = 500  mUligauss  and  for  PFRR  Bq  = 560  milligauss).  Plots  of  pitch 
angle,  defined  as  180°  - 7,  for  each  rocket  as  a function  of  rocket  altitude  are 
shown  in  Figures  18,  19,  20  and  21.  These  data  represent  the  basic  information 
from  which  rocket  aspect  relative  to  local  vertical  (zenith  angle)  is  deduced  in  the 
discussion  section  of  this  report.  The  magnetometer  data  from  rockets  A30. 205-5 
and  6 was  multiplexed  with  horizon  sensor  data  causing  degraded  magnetometer 
records  that  could  not  be  processed  with  standard  computer  techniques.  The  data 
shown  in  Figures  18  and  13  for  these  rockets  were  manually  reduced. 

In  considering  the  magnetometer  data  presented,  it  is  important  to  understand 
that  the  pitch  angle  as  computed  represents  a cone  about  the  magnetic  field  vector 
on  which  the  rocket  axis  can  lie.  Any  position  on  this  cone  will  produce  the  same 
magnetometer  output  voltage.  For  this  reason  it  is  essential  to  introduce  addi- 
tional evidence  and/or  assumptions  to  deduce  the  rocket  position  on  the  locus  cone. 
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5.  DISCUSSION 


5.1  Intensity  Ratio  Profiles 

From  an  inspection  of  the  measured  emission  profiles  shown  in  Figures  18, 

19,  20  and  21,  it  is  clearly  evident  that  the  data  are  significantly  modulated  by  the 

rocket  coning  motion.  In  general,  this  can  be  attributed  to  a van  Rhijn  effect 
8 

(Peterson  ) if  the  emissions  are  assumed  to  emanate  from  a homogeneous  ootically 
thin  layer.  The  effects  of  aspect  must  be  accounted  f^r  if  absolute  values  of 
intensity  are  to  be  compared  between  rockets.  This  is  usually  accomplished  by 
correcting  all  data  to  zenith  profiles  using  a van  Rhijn  factor,  viz. 


R 


o 


(5) 


where 

K = zenith  radiance  at  altitude  h 
o 

Rq  = radiance  measured  at  an  angle  a from  the  vertical  at  an  altitude  h 

Vq  = van  Rhijn  factor  for  angle  a at  an  altitude  h. 

For  the  two  channel  radiometer  systems  used  in  these  measurements,  it  can  be 
assumed  that  the  van  Rhijn  factors  are  the  same  for  each  channel  on  any  one  rocket 
because  both  have  the  same  aspect  at  any  instant  in  time  and  the  measured  emis- 
sions emanate  from  about  the  same  altitude  (Rogers  et  al*).  Thus,  the  ratio  of 
measured  intensities  from  the  two  channels  on  any  one  rocket  is  the  same  as  the 
ratio  of  the  zenith  emissions.  Consequently,  the  ratios  of  the  two  wavelength 
channels  from  each  rocket  measurement  can  be  compared  with  the  ratios  measured 
with  the  other  rockets.  The  ratio  of  the  short  wavelength  channel  (1.  5642  pm)  to 
the  long  wavelength  channel  (1.  985  nm)  on  each  rocket  were  computed  and  plotted 
as  a function  of  altitude  (see  Figure  22). 

Considering  first  the  ratio  data  shown  in  Figure  22  for  emission  intensities 
measured  on  rocket  AO30. 311-3  (midlatitude  nighttime  flight),  there  is  clearly 
little  effect  of  aspect  compared  to  the  aspect  modulation  on  the  individual  channel 
intensity  data  in  Figure  21.  The  data  for  this  rocket  show  a nearly  constant  ratio 
value  of  about  1.4  up  to  83  km.  As  the  rocket  penetrates  the  main  layer  of  emission 
(above  83  km),  the  ratio  diminishes  to  a value  of  about  0.  5 at  93  km.  At  92.  5 km, 
the  measured  long  and  short  wavelength  intensities  are  equal  as  the  ratio  is  1.  The 


8.  Peterson,  A.W.,  and  Kieffaber,  L.M.  (1972)  A irglow  fluctuations  at  2. 2u, 
J.  Atmos.  Terr.  Phys.  34:1357. 
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data  shown  for  rockets  A30,  205-5  and  A 30.  205-6  (high  latitude  nighttime  flights) 
are  in  agreement,  showing  a slight  increase  with  altitude  in  the  ratio,  with  values 
from  l.  2 to  1.  6.  These  results  do  not  agree  with  ratio  data  reported  by  Rogers 
et  al*  for  high  latitude  nighttime  measurements  where  the  long  wavelength  meas- 
ured intensities  were  higher  than  the  short  wavelength.  The  discrepancy  could  be 
a result  of  slightly  different  spectral  response  characteristics  of  the  Rogers'  meas- 
urement, even  though  the  reported  spectral  passbands  (1.4  to  1.  65  pm  and  1.  85  to 
2.  12  pm)  seem  to  be  the  same  as  those  reported  here.  The  discrepancy  is  more 
likely  due  to  the  uncertainties  in  that  first  measurement,  which  was  really  a flight 
test  of  the  instrument  in  which  lack  of  time  precluded  a precise  instrument  calibra- 
tion. 

The  ratio  data  in  Figure  22  for  rocket  AO30.  311-2  (midlatitude  twilight  flight) 
show  different  characteristics  than  the  nighttime  data.  These  data  show  a higher 
ratio  and  considerably  more  variability.  Since  the  rocket  was  flown  at  a time  when 
the  earth  shadow  height  was  about  60  km  (x  = 99°),  data  above  this  altitude  were 
probably  periodically  contaminated  by  sunlight  reflecting  into  the  instrument 
depending  upon  the  vehicle  aspect.  This  seems  to  be  confirmed  by  the  fact  that  the 
short  wavelength  intensity  profile  (Figure  20)  shows  a greater  modulation  with 
aspect  than  the  long  wavelength  profile.  In  any  case,  ratio  data  below  60  km  should 
be  unaffected  by  sunlight.  In  this  altitude  region,  the  ratio  is  3.4.  If  it  is  now 
assumed  that  the  minimum  ratio  values  above  60  krn  are  valid,  the  ratio  profile 
takes  on  a characteristic  of  diminishing  values  with  altitude  (dashed  line  in  Figure 
22)  thus  implying  a markedly  different  OH  production  mechanism  during  twilight 
than  at  nighttime. 

5.2  Volume  Kmission  Rates 

In  order  to  derive  volume  emission  rates  from  the  rocket  measured  data,  it 
is  necessary  to  correct  the  measured  altitude  profiles  for  aspect  (compute  zenith 
profiles)  and  then  differentiate  these  profiles  with  respec;  to  altitude,  viz. 


(6) 


where 

N = volume  emission  rate  (megaphotons/ cc-sec) 
Rq  = zenith  radiance  (rayleighs) 

h = altitude  (cm). 
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This  method  is  applicable  only  for  profiles  that  traverse  the  emitting  region  and 
for  emission  regions  which  are  homogeneous,  two-dimensionally  optically  thin, 
and  void  of  extinction  effects  such  as  absorption  (van  Rhijn  regions).  The  emission 
criteria  can  be  assumed  satisfied  for  all  the  hydroxyl  data  presented,  but  only 
rocket  AO30.  311-3  provided  measurements  through  the  principle  hydroxyl  emis- 
sion region.  Consequently,  only  the  two  emission  profiles  measured  from  this 
rocket  were  considered  for  computing  volume  emission  rates. 

The  first  step  in  the  data  processing  is  to  derive  zenith  profiles  from  the 
measured  data.  This  requires  computation  of  the  instrument  viewing  angles  (p) 
relative  to  the  local  vertical  for  each  measured  radiance  data  point.  Unfortunately, 
the  only  direct  rocket  aspect  data  available  is  that  from  the  magnetometer  and  these 
data  alone  are  not  sufficient  to  derive  the  aspect  angle  a.  However,  the  measured 
radiance  data  also  contains  some  aspect  information.  For  instance,  modulation  of 
the  measured  radiance  data  in  Figure  21  below  the  main  emission  region  is  indica- 
tive of  changes  in  rocket  aspect.  In  fact,  the  modulation  of  the  radiance  data  show 
an  in-phase  correlation  with  the  magnetometer  pitch  angle  data.  If  one  considers 
an  idealized  case  where  a rocketborne  radiometer  is  flown  vertically  (with  no  coning) 
up  to  and  through  a van  Rhijn  type  emission  layer,  the  measured  radiance  would  be 
a zenith  profile.  The  measured  radiance  level  would  be  constant  up  to  the  bottom 
of  the  region  and  then  would  diminish  as  the  instrument  traversed  the  layer.  Thus, 
for  an  actual  measurement  case  where  coning  and  non-vertical  trajectory  exist  and 
it  can  be  assumed  that  there  are  no  or  at  least  negligible  hydroxyl  emissions  below 
some  altitude,  a measured  emission  profile  that  has  been  properly  corrected  for 
aspect  (zenith  profile)  should  show  a constant  radiance  level  below  the  layer.  A 
technique  based  on  this  idea  was  devised  to  correct  the  measured  profiles  from 
rocket  AO30.  311-3.  In  this  approach,  the  magnetometer  data  were  used  to  compute 
viewing  angles  (a)  by  assuming  various  headings  for  the  rocket  axis  of  precession. 
These  aspect  angles  were  then  used  to  correct  the  measured  profiles  for  each  head- 
ing using  the  van  Rhijn  method,  The  corrected  profiles  were  plotted  and  compared 
for  the  "best  fit"  to  constant  radiance  below  75  km,  and  this  best  fitted  profile  was 
then  differentiated  by  a Fourier  technique  to  obtain  volume  emission  rates.  The 
details  of  this  technique  follow. 

The  original  measured  radiance  data  files  (1660  data  points  per  second)  were 
averaged  to  obtain  one  data  point  every  0.  5 km,  correlated  with  altitude  and  pitch 
angle.  These  data  are  plotted  in  Figure  23.  The  pitch  angle  (P.  A. ) data  were 
computed  using  the  magnetometer  output  and  Eq.  (4)  to  find  angle  y and  then  pitch 
angle  from 

P.  A.  = 180°  - y . (7) 
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Figure  23.  Smooth  Radiance  Data  for  the  1.  5642-nm  and  the  1.  985-pm 
Channels  Measured  on  Rocket  AO30  311-3.  Data  not  corrected  for  rocket 
aspect 


The  angle  (9^)  between  the  rocket  axis  of  precession  and  the  earth's  magnetic  field 
was  determined  from  the  envelope  of  y as  follows  (see  Appendix  Figure  C2): 


e 

c 


y + y . 

_ 'max  'mm 

2 


(8) 


The  rocket  axis  of  precession  is  defined  in  space  by  angle  (?c  if  the  heading  (9  ) of 
the  axis  is  known  (the  position  of  the  earth's  magnetic  field  is  assumed  to  be  known). 
The  angle  «c  between  the  local  vertical  and  the  axis  of  precession  was  computed 
for  several  choices  of  ¥c  and  each  data  point  from  (see  Appendix  C) 
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where 

1 = cos  9 

c 

m = cos  6 cos  (D  - 
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n = sin  6 

V = heading  relative  to  true  north  of  axis  of  precession 
6 = earth  magnetic  field  dip  angle  (60°  for  WSMR) 

D = earth  magnetic  field  declination  (+10°  for  WSMR). 

The  principal  assumption  in  this  calculation  is  that  the  heading  of  the  rocket  axis 
of  precession  does  not  change  during  the  measurement  period  from  about  53  km  to 
106  km.  Since  there  is  little  air  drag  in  this  region  and  no  other  perturbations 
such  as  tip  ejections  taking  place,  the  assumption  is  presumed  to  be  valid. 

The  aspect  angle  (p ) between  the  local  vertical  and  the  rocket  axi3  is  derived 
in  detail  in  Appendix  C and  is  given  by 

a * oc  + - y . (10) 

Equation  (10)  is  an  approximation  based  on  the  assumptions  that  the  radiance  data 
modulation  due  to  coning  is  in-phase  with  pitch  angle  and  that  the  half  cone  angle 
(0)  is  small  such  that  ccs  0 ~ 1 and  sin  0 * 0.  Both  of  these  assumptions  were 
established  appropriate  for  rocket  AO30.  311-3. 

The  van  Rhijn  correction  factors  (V)  were  computed  from  Eq.  (11)  for  each 
data  point  and  applied  to  the  measured  data  with  Eq.  (5)  where 


V * van  Rhijn  factor  for  a thin  layer  at  z ^ 
V2  = van  Rhijn  factor  for  a thin  layer  at  x2 
R^.  * radius  of  earth 
h - height  above  earth  of  measurement 
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Zj  = altitude  of  lower  !■  id  of  emitting  layer  = 76  km 

z^  « altitude  of  upper  bound  of  emitting  layer  = 92  kin. 

In  the  computation  of  V the  rocket  altitude  (h)  replaces  the  quantity  z^  when 
h > z j,  so  that  V j = sec  a whenever  the  rocket  is  within  the  emitting  layer.  Several 
zenith  profiles  were  thus  computed,  one  for  each  assumed  heading  ('I'  ).  The  zenith 
profile  having  the  best  fit  to  constant  radiance  below  75  km  is  shown  in  Figure  24. 
Also  shown  is  aspect  angle  a. 
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Figure  24.  Zenith  Altitude  Radiance  Profiles  of  the  1.  5642 -pm  and  the 
1.  985-pm  Channels  Computed  from  Data  Measured  on  Rocket  AO30.  311-3. 
Also  shown  is  rocket  aspect  angle  (p)  relative  to  the  local  vertical 


The  modulation  still  present  in  the  zenith  profile  below  the  emission  layer 
after  the  van  Rhijn  correction  (see  Figure  24)  can  be  attributed  to  two  assumptions 
used  in  deriving  Eq.  (10).  First  it  is  assumed  that  the  half  cone  angle  (0)  is  con- 
stant throughout  the  flight,  and  second  that  the  modulation  of  the  uncorrected  data 
is  inphase  with  the  pitch  angle  (P.  A. ).  The  error  in  the  aspect  angle  <p)  when 
Eq.  (10)  is  used,  is  dependent  on  the  half  cone  angle  (0)  and  on  the  difference  in 
heading  between  the  rocket  axis  (&c)  and  the  heading  of  the  local  earth  magnetic 
field  (D).  The  relationship  for  the  maximum  error  (emaJ£)  is  derived  in  Appendix  C 
and  is  given  by 
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and 

e * maximum  error  in  a (degrees) 

max 

ji  * half  cone  angle  (less  than  5°  below  76  km  and  6°  above  76  km) 

Q * phase  angle 

1t,c  * heading  of  the  axis  of  precession  (300°) 

D * earth  magnetic  field  declination  (10°  at  WSMR) 

6 ■ earth  magnetic  dip  angle  (60°  at  WSMR) 

0c  ■ angle  between  axis  of  precession  and  earth  magnetic  field. 

Below  the  layer,  6 ■ 143°  (see  Eq.  (8))  so  that  emax  * ±4.  3°.  This  means  that  the 
angle  a used  to  correct  the  profile  was  too  large  by  4.  3°  at  the  minimum  radiance 
values  in  Figure  24  and  too  small  by  4.  3°  for  the  peak  radiances  shown  below 
76  km.  These  corrections  were  applied  to  the  peak  and  minimum  data  points  of  the 
radiance  modulations  in  Figure  24.  The  resulting  data  points  fell  within  1 KR  of  the 
mean  values  of  148. 4 KR  for  the  1.  56-Mrn  channel  and  104. 4 KR  for  the  1.  99-pm 
channel.  Consequently,  the  axis  of  precession  heading  of  300°  was  concluded  to  be 
correct.  The  same  technique  was  applied  to  the  radiance  data  within  the  layer  where 

8 * 131°  and  B was  nearly  constant  at  6°.  In  this  case,  e * ±3.  95°  and  the 

radiance  values  corresponding  to  the  peak  and  minimum  pitch  angles  were  similarly 
corrected.  It  was  found  that  the  actual  radiance  errors  were  only  of  the  order  of  a 
few  kilorayleighs.  Thus,  it  was  concluded  that  the  radiance  data  within  the  layer 
were  adequately  corrected  by  the  initial  approximation  of  a computed  from  Eq.  (10). 
The  final  zenith  emission  profiles  were  deduced  by  this  analyses  and  are  shown  in 
Figure  24  by  the  solid  lines  where  the  error  bar  below  76  km  is  9 KR,  and  above 
76  km  the  error  bar  is  about  4 KR. 

Computation  of  volume  emission  rates  by  differentiating  zenith  profile  data  as 
defined  by  Eq.  (6)  requires  sophisticated  techniques  when  digital  data  are  involved. 
The  difficulty  arises  from  the  fact  that  normal  random  fluctuations  between  consec- 
utive data  points  due  to  teleme«,ry  and  recorder  noise,  give  rise  to  large  slope 
changes  when  simple  differentiation  is  employed.  Fitting  the  data  with  continuous 
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functions  such  as  polynomials  provide  the  necessary  smoothing  of  the  data  points, 
but  this  method  can  overly  influence  the  derivative  and  thus  ihe  volume  emission 

9 

rate  computed,  depending  on  the  function  used.  Another  technique  used  by  Pearce 
and  also  by  Rogers*  is  to  smooth  the  data  by  Gaussian  weights  piior  to  differentia- 
ting. This  method,  though  successful,  does  not  take  advantage  of  discriminating 
against  data  fluctuations  beyond  the  frequency  response  of  the  measuring  instrument. 
The  data  processing  technique  selected  for  results  reported  here  employs  deriva- 
tive digital  filtering  and  Fast  Fourier  transform  to  delete  data  fluctuations  outside 
the  instrument  frequency  response  (about  1 Hz)  and  compute  the  derivative  dRQ/dt. 
The  volume  emissions  rate  (N)  was  then  computed  by  the  chain  rule 
N = -(dR  /dt)(dt/dh).  A complete  description  of  the  method  is  given  in  Appendix  D. 
In  order  to  provide  a large  data  base  for  the  derivative  digital  filtering  process,  the 
initial  data  base  consisting  of  1660  points  per  second  was  used.  Aspect  correc- 
tions were  applied  using  the  van  Rhijn  Eqs.  (5)  and  (11),  and  the  viewing  ange  (») 
derived  previously. 

The  volume  emission  rates  were  derived  in  this  manner  for  the  two  radiometer 
channels  on  rocket  AO30.  311-3  and  are  shown  in  Figure  25.  The  1.  56-/nm  emission 


Figure  25.  Volume  Emission  Rates  for  the  1.  5645-um 
and  1.985-um  Channels  Computed  From  Data  Measured 
on  Rocket  AO30.  311-3 


9.  Pearce,  J.  B.  (1969)  Rocket  measurements  of  nitric  oxide  between  60  and  95 
kilorr  .“ters,  J.  Geophys.  Res.  74:853. 
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5 -3-1 

rate  shows  a peak  of  1.  9 X 10  photons  cm  sec  at  89  km,  while  the  1.  99-pm 

5 -3  -1 

emission  rate  peaks  at  1,  5 x io  photons  cm  sec  at  89.  5 km.  Most  of  the 

emission  in  both  wavelength  bands  come  from  a narrow  layer  8-  to  10-km  thick, 

with  the  largest  difference  in  volume  emission  rates  occurring  at  and  below  the 

altitude  of  peak  emission.  In  this  region,  the  short  wavelength  emissions  exceed 

4 -3-1 

the  long  wavelength  emissions  by  approximately  6 X 10  photons  cm  sec  . 
Above  the  peak  of  the  layer,  the  two  emissions  are  nearly  equal. 

If  the  instrument  spectral  bandpasses  were  flat  and  the  emissions  were  con- 
stant throughout  these  spectral  regions,  the  peak  volume  emission  rates  (N)  per 
micrometer  for  each  channel  would  be  approximately 

-3  - 1 - 1 

N(l.  56  pm)  = 1 megaphoton  cm  sec  pm 

-3  - 1 - 1 

N(l,  99  pm)  = 0.  7 megaphoton  cm  sec  pm 


6.  CONCLUSIONS 

The  following  conclusions  are  based  solely  on  the  rocketborne  radiometer 
measurements  made  in  the  1.  56-pm  (AX  = 0.  2 1 pm)  and  the  1.  99-pm  (AX  = 0.  20  pm) 
infrared  bands  as  described  in  this  report: 

(1)  The  principle  hydroxyl  emissions  in  the  bands  stated  occur  in  a thin  layer 
8-  to  10-km  thick  with  peak  volume  emission  rates  occurring  at  approximately 

89  km. 

(2)  There  is  no  conclusive  evidence  to  substantiate  any  appreciable  OH  emis- 
sion layers  below  80  km  during  night  conditions  in  the  bands  measured, 

(3)  The  ratio  of  the  nighttime  1.  56  pm  to  the  1.  99  pm  overhead  radiance  is 
constant  at  a value  of  approximately  1.4  up  to  80  km  regardless  of  latitude. 

(4)  There  is  some  evidence  that  at  twilight  the  ratio  of  th  • 1.  56  pm  to  the 
1.  99  pm  radiance  is  larger  than  1.  4 and  diminishes  with  altitu  K up  to  at  least 
80  km. 

(5)  Estimated  nighttime  peak  volume  emission  rates  (N)  based  on  one  rocket 
flight  at  WSMR  Pie- 

-3  - 1 - 1 

N(l.  56  pm)  = 1 megaphoton  cm  sec  pm 

-3  - 1 - 1 

N(l.  56  pm)  = 0.7  megaphoton  cm  sec  pm 

An  attempt  has  been  made  in  this  report  io  document  in  detail  the  techniques 
used  in  measuring  and  reducing  the  rocket  radiometer  data.  On  the  other  hand, 
little  emphasis  has  been  placed  on  hydroxyl  chemistry  or  comparing  the  reduced 
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hydroxyl  radiance  and  volume  emission  rate  data  with  those  reported  by  other 
workers  in  the  field.  The  reasons  for  this  approach  were: 

(1)  Hydroxyl  radiometric  data  reported  to  date  were  measured  in  spectral 
bandpasses  sufficiently  different  from  those  reported  here  such  that  direct  com- 
parison has  little  meaning. 

(2)  In  much  of  the  radiance  data  reported  in  the  literature,  there  is  insufficient 
information  on  the  measuring  instrument  spectral  responses  to  accurately  assess 
differences  in  reported  results. 

(3)  A proper  analysis  of  chemical  processes  was  considered  beyond  the  scope 
of  this  report,  because  it  involves  comparison  of  the  measured  data  with  radiances 
found  by  integration  of  known  or  modeled  spectral  OH  emissions  throughout  the 
spectral  bandpasses  of  the  measured  data  as  well  as  consideration  of  reaction  rates 
and  upper  atmospheric  composition.  The  authors  feel  that  sufficient  information 
has  been  presented  in  this  report  to  facilitate  the  accomplishment  of  the  chemistry 
study  mentioned  above,  but  this  should  be  the  subject  of  another  report. 
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Appendix  A 

Selected  Rocket  Instrument  Parameters 


I.  OH  RADIOMETER  FIGURES  OF  MERIT 

This  appendix  includes  selected  figures  of  merit  describing  the  OH  radiometers 
(Tables  A 1 through  A4)  used  on  Astrobee  D rockets  A30.  205-5,  A30.  205-6, 

AO30.  311-1,  AO30.  311-2,  and  AO30.  311-3  (Figures  A 1 and  A2  for  AO30.  311-3). 
The  primary  reference  sources  for  these  data  are  a series  of  calibration  manuals 
provided  by  Wyatt.  The  data  presented  applies  to  all  the  rocket  instruments  in 
general  except  where  noted.  For  example,  the  instrument  relative  spectral  re- 
sponse is  the  same  for  all  rocket  instruments,  whereas  fields-of-view  are  different 
for  each  rocket. 
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Table  A 1 


OH  RADIOMETER 

Channel  2 

1.985  MICRONS 

THE  INSTRUMENTS 

RELATIVE  SPECTRAL 

RESPONSE  CURVE 

WAVELENGTH 

RELATIVE 

WAVELENGTH 

RELATIVE 

(MICRONS) 

RESPONSE 

(MICRONS) 

RESPONSE 

1.7600 

.00 

1.7700 

.00 

1.7800 

.01 

1.7900 

.01 

1.8000 

.01 

1.3100 

.04 

1.8200 

.10 

1.8300 

.13 

1.8400 

.19 

1.8500 

.40 

1.8600 

.60 

1.8700 

.46 

1.8800 

.46 

1.8900 

.41 

1.9000 

.71 

1.9100 

.75 

1.9200 

.72 

1.9300 

.75 

1.9400 

.76 

1.9500 

.68 

1.9600 

.53 

1.9700 

.50 

1.9800 

.56 

1.9900 

.77 

2.0000 

.66 

2.0100 

1.00 

2.0200 

.84 

2.0300 

.72 

2.0400 

.73 

2.0500 

.83 

2.0600 

.83 

2.0700 

.69 

2.0800 

.63 

2.0900 

.70 

2.1000 

.81 

2.11C0 

.80 

2.1200 

.58 

2.1300 

.26 

2.1400 

.08 

2.1500 

.02 

2.1600 

.01 

2.1700 

.01 

2.1800 

.01 

2.1900 

.01 

2.2000 

.01 

2.2100 

.01 

2.2200 

.01 

2.2300 

.00 

Table  A2 


OH  Radiometer 

Channel  1 

1.5642 

MICRONS 

THE  INSTRUMENTS  RELATIVE 

SPECTRAL  RESPONSE 

CURVE 

WAVELENGTH 

RELATIVE 

WAVELENGTH 

RELATIVE 

(MICRONS) 

RESPONSE 

(MICRONS) 

RESPONSE 

1.3000 

.00 

1.3100 

.00 

1.3200 

.00 

1.3300 

.01 

1.3400 

.01 

1.3500 

.01 

1.3600 

.02 

1.3700 

.02 

1.3800 

.03 

1.3900 

.04 

1.4000 

.06 

1.4100 

.09 

1.4200 

.13 

1.4300 

.20 

1.4400 

.29 

1.4500 

.39 

1.4600 

.47 

1.4700 

.67 

1.4800 

.84 

1.4900 

.90 

1.5000 

.92 

1.5100 

.95 

1.5200 

.96 

1.5300 

.92 

1.5400 

.84 

1.5500 

.84 

1.5600 

.85 

1.5700 

.87 

1.5800 

.91 

1.5900 

.93 

1.6000 

.93 

1.6100 

.95 

1.6200 

.90 

1.6300 

.32 

1.6400 

.84 

1.6500 

.92 

1.6600 

1.00 

1.6700 

.47 

1.6800 

.11 

1.6900 

.02 

1.7000 

.01 

1.7100 

.01 

1.7200 

.01 

1.7300 

.00 
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Table  A3 


NR-3B-11  CHANNEL  1 

1.985  MICRON 

1.985  MICRON 

FILTER  RESPONSE 

WAVELENGTH 

RELATIVE 

WAVELENGTH 

RELATIVE 

(MICRONS) 

RESPONSE 

(MICRONS) 

RESPONSE 

1.7600 

.00 

1.7700 

.40 

1.7800 

.50 

1.7900 

.70 

1.8000 

1.20 

1.8100 

4.00 

1.8200 

9.70 

1.8300 

12.00 

1.8400 

17.00 

1.8500 

37.00 

1.8600 

54.80 

1.8700 

41.50 

1.8800 

41.00 

1.8900 

36.90 

1.9000 

63.00 

1.9100 

67.00 

1.9200 

63.50 

1.9300 

66.00 

1.9400 

66.20 

1.9500 

58.50 

1.9600 

45.50 

1.9700 

42.50 

1.9800 

47.20 

1.9900 

61.00 

2.0000 

80.00 

2.0100 

83.50 

2.0200 

70.00 

2.0300 

59.50 

2.0400 

60.00 

2.0500 

68.00 

2.0600 

67.50 

2.0700 

56.00 

2.0800 

51.00 

2.0900 

56.00 

2.1000 

64.00 

2.1100 

62.40 

2.1200 

45.00 

2.1300 

20.00 

2.1400 

6.00 

2.1500 

1.70 

2.1600 

1.00 

2.1700 

.70 

2.1800 

.70 

2.1900 

.70 

2.2000 

.60 

2.2100 

.60 

2.2200 

.40 

2.2300 

.00 

Table  A4 


NR-3B-1X  CHANNEL  2 

1.5642  MICRON 

1.5624  MICRON 

FILTER  RESPONSE 

WAVELENGTH 

RELATIVE 

WAVELENGTH 

RELATIVE 

(MICRONS) 

RESPONSE 

(MICRONS) 

RESPONSE 

3.6780 

.00 

3.6880 

.02 

3.6980 

.05 

3.7080 

.08 

3.7180 

.10 

3.7280 

.13 

3.7380 

.20 

3.7480 

.24 

3.7580 

.35 

3.7680 

.43 

3.7780 

.60 

3.7880 

.90 

3.7980 

1.30 

3.8080 

1.96 

3.8180 

2.80 

3.8280 

3.76 

3.8380 

4.50 

3.8480 

6.40 

3.8580 

8.00 

3.8680 

8.62 

3.8780 

8.79 

3.8880 

9.00 

3.8980 

9.08 

3.9080 

8.70 

3.9180 

8.00 

3.9280 

7.87 

3.9380 

8.00 

3.9480 

8.10 

3.9580 

8.42 

3.9680 

8.54 

3.9780 

8.47 

3.9880 

8.60 

3.9980 

8.05 

4.0080 

7.30 

4.0180 

7.40 

4.0280 

8.02 

4.0380 

8.65 

4.0480 

4.00 

4.0580 

.90 

4.0680 

.20 

4.0780 

.10 

4.0880 

.06 

4.0980 

.05 

4.1080 

.00 
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Appendix  B 

Tabulations  of  Rodent  AO 30. 31 1-3  Data 


This  appendix  contains  tabulations  of  radiance  and  rocket  aspect  data  as  a 
function  of  rocket  altitude  for  rocket  AO30.  311-3.  The  table  labeled  "Averaged 
Data"  consists  of  an  average  of  all  data  points  approximately  every  0.  5 km  with 
spurious  data  points  removed.  These  data  are  the  input  used  in  computation  of  the 
data  shown  in  the  table  labeled  "Zenith  Profile.  " Symbol  definitions  are  as  follows: 

ALTITUDE  - Rocket  altitude  (kilometers) 

PITCH  - Rocket  pitch  angle  (P.  A. ) (degrees) 

IRIG  16  - Radiance  in  the  1.  56  pm  radiometer  channel  (kilorayleighs) 

IRIG  15  - Radiance  in  the  1.  99  pm  radiometer  channel  (kilorayleighs) 

RATIO  - Ratio  of  radiance  IRIG  16  to  IRIG  15 

ALPH  - Rocket  aspect  (viewing  angle  of  radiometer)  relative  to  the  local 
vertical  (degrees) 
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AO  30.311-1  ASCENT 
AVERAGEO  flATA 
IRIG  16 
TRIG  15 


ALTITUOE 

56.966 

57.446 
56*020 

58.497 
56.970 
59.53*. 
60.  00  3 

60. 469 
61.0  2'. 
61.485 
61.943 
62.490 
63.034 

63.483 
63.931 
64.466 

64.996 

65.523 
66.049 
66.569 
66.999 
67.426 
67.939 

68.447 

68.95? 

69.45* 

69.950 

70.44*= 

70.935 

71.50* 

72.069 

72.547 

73.023 

73.496 

73.966 
74.431 
75.656 
75.95° 

76.484 
77.741 
78.03' 

78.469 
78.97? 
79.471 
79.964 

80.523 
81.007 

81.485 
82.029 

82.497 


PITCH 

35.416 
33.501 
32.551 
32.729 

34.772 
37. 758 
40.091 
41.702 
42.634 
41.734 
40.697 
38.543 
36.435 
35.474 
34.548 
36.521 
34.875 

35.417 
36.178 
36. 667 
37.101 
37.188 
37.346 
17.146 
36.455 
35.907 
35.131 
34.  894 
34.762 
35.663 
37.218 
39.040 
41 . 358 

43.772 
45.793 
4C.004 
51.816 
52.916 
53.336 
52.353 
51.910 
50.853 
49.591 
48.085 
46.636 
45.030 
44.H2 
43.725 
43.639 
44.055 


IRIG16 

151.503 

150.144 

148.260 
147.535 
145.940 
150.289 
154.112 
156.559 
160.327 
15°. 4°4 
158.443 
155.254 
153. 224 
148. 622 
148.803 
148. 167 
146.919 
147.426 

148.948 
150.  108 
15  0.796 

153.261 
154.239 
158.226 
158.044 
156.776 
157.501 
155.870 
155.109 
155.281 
156.124 
156. 743 
159.349 
160.192 
161.916 
165.256 
176.790 
179.643 
182.498 
183.412 
191.636 

180.948 
176.046 
173.881 
17C.937 
167.  ->50 
163. 118 
161.668 
15°.902 

152. 246 


IRIG15 

112.070 

103.853 

105.497 

104.246 
103.639 

108.497 
108.104 

110.641 

111.641 
113.999 

109.140 
106.068 
107.712 

104.139 

104.854 

102.996 

105.711 

104.997 
107.569 

106.711 

107.247 

107.712 

106.140 
109.784 
109.355 
109.355 
107.925 

107.497 

199.140 

107.783 

106.783 

109.283 

107.283 
112.856 

114.285 

117.285 
135.035 
126.217 
128.199 
126.645 
127.717 
125.002 
123.216 
120.501 
118.768 
116.142 
113. 9^1 
113. 1-.? 
110.730 
108.648 
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82.96? 

44.818 

151. 340 

108. 319 

83.48^ 

46.365 

147.236 

109.391 

84.007 

48.193 

143. 213 

109.498 

84.51Q 

49.792 

140.871 

108.640 

85.027 

51.607 

133.891 

105.747 

85.527 

52.620 

127.150 

107.622 

86 . ( 2? 

53.763 

122. 258 

105.640 

86.5 il 

54.492 

115.408 

102.639 

86 . f 94 

55.032 

107. 554 

100.281 

87.470 

54.935 

98.231 

95.458 

87.999 

54.315 

89.044 

92.865 

88 .521 

53.484 

76. 378 

84. 741 

88.977 

52. 44  7 

66.458 

76.382 

89.484 

50. 842 

53.123 

67.615 

89.983 

48.843 

41.585 

56.876 

90.474 

46. 855 

34.707 

49.696 

91.011 

44.946 

28.944 

42.365 

91.539 

43. 368 

24. 443 

36.235 

92. 006 

42. 640 

21. 557 

32.870 

92.464 

42.391 

15. 2<*6 

31.177 

92.965 

43.057 

17.681 

28.004 

93.505 

44.583 

17.181 

24.796 

94.033 

46.652 

15.898 

25.132 

94.504 

43  « 868 

15.702 

22.774 

95.550 

53.014 

13.349 

20.277 

95.988 

54. 153 

15.025 

21.009 

96.503 

54.793 

12.615 

20.581 

97.003 

54. 640 

12.579 

15.758 

97.490 

53.658 

12. 252 

15.9J7 

98.002 

52.111 

12. 366 

14.513 

98.498 

49. 773 

12.  180 

14.365 

98.978 

46.888 

12.180 

14.911 

99.477 

44.047 

12.180 

13.196 

99.857 

43.140 

12.180 

12.400 

ii 

I 
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AO  30.311-'?  ASCENT 
ZENITH  PROFTLES 
IRIG16 
IRIG15 

HEADING  = TOO.  OEG. 


ALTITUOE 

IPIG16 

56.966 

147.934 

57.446 

147.608 

58.020 

146.185 

58.497 

145. 394 

58.970 

142.847 

59.634 

145. 301 

60.003 

147.271 

60.469 

148. 252 

61.024 

150.963 

61.486 

151.004 

61.943 

150.907 

62.490 

149.544 

63.031* 

149.004 

63.48? 

145. 089 

63.931 

145. 767 

64.466 

145.178 

64.996 

143.751 

65.52? 

143.953 

66.048 

144.999 

66.568 

145.892 

66.999 

146.224 

67.428 

148.556 

67.939 

149.400 

68.447 

153. 397 

68.952 

153.677 

69.453 

152.787 

69.950 

153.958 

70.446 

152. 369 

70.935 

151.484 

71.504 

150. 856 

72,068 

150. 344 

72.547 

149.269 

73.023 

149.391 

73.496 

147. 397 

73.965 

145.542 

74.431 

14: . 208 

75.656 

149.602 

75.959 

149.916 

76.484 

151.131 

77.741 

152.705 

78.03? 

151.797 

78.469 

152.741 

78.97? 

150.342 

79.471 

150.722 

79.964 

150.277 

80.523 

149. 034 

81.007 

146.688 

81,486 

145. 858 

82.028 

144.368 

82.497 

136.978 

82.96? 

135. 270 

IRIG15 

RATIO 

109.430 

1.352 

102.099 

1.446 

104.021 

1.405 

102.733 

1.415 

101. 443 

1.408 

104.897 

1.385 

103.306 

1.426 

104.771 

1.415 

105. 121 

1.436 

107.930 

1.399 

103.950 

1.452 

102. 167 

1.464 

10  4,74*. 

1.423 

101.663 

1.427 

102.714 

1.419 

100.905 

1.439 

103.432 

1.390 

102.523 

1.404 

104.717 

1.385 

103.  714 

1.407 

103.995 

1.406 

104.405 

1.423 

102.809 

1.453 

106,433 

1.441 

106.333 

1.445 

106.573 

1.434 

105.498 

1.459 

105. 083 

1.450 

106.590 

1.421 

104,711 

1.441 

102. 830 
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Appendix  C 

Racket  Aspect  From  Magnetometer  Data 


A magnetometer  responds  to  the  effective  magnetic  field  component  parallel  to 
its  sensing  axis  as  follows: 


v 


= y + JL 
Vb  k 


(Cl) 


where 

v = voltage  output  from  magnetometer  due  to  field  B 

3 = instantaneous  effective  magnetic  field  intensity  component  aligned  with 
the  magnetometer  axis  (milligauss) 

= bias  voltage  (vclts)  (voltage  output  from  magnetometer  when  B = 0) 

k = calibration  factor  (milligauss/volt). 

Magnetometers  are  calibrated  by  applying  known  fields  to  the  instruments  in  a 
controlled  field  environment  and  measuring  the  output  voltage.  The  bias  voltage 
and  calibration  factor  k are  determined  from  these  data  (see  Figure  Cl).  Equa- 
tion (Cl)  is  valid  only  in  the  linear  portion  of  the  calibration  curve  shown  in 
Figure  Cl. 

In  rocket  research,  magnetometers  are  used  in  the  payloads  to  determine 
rocket  spin  and  precession  (coning)  relative  to  the  earth's  magnetic  field.  In  this 
application,  the  contribution  to  the  effective  magnetic  field  acting  on  the  magneto- 
meter from  the  earth's  magnetic  field  is 
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Figure  Cl.  Typical  Magnetometer  Calibration  Data 
(Schonstedt  Instrument  Company) 


B2  = Bq  cos  <t>  (C2) 

where 

Bj  = component  of  earth's  field  aligned  with  the  magnetometer  axis 
Bq  = total  field  intensity  of  earth's  magnetic  field 

$ = instantaneous  angle  between  earth's  total  field  and  magnetometer  axis 

Let  AB  be  the  field  intensity  aligned  with  the  magnetometer  axis  arising  from  cur- 
rents in  the  payload  so  that  the  net  effective  field  B acting  on  the  magnetometer  is 

B = Bl  + AB  = Bq  cos  <t>  + AB  . (C3) 

Combining  Eq.  (C3)  and  Cl)  yields  the  instantaneous  voltage  output  from  the  mag- 
netometer due  to  all  fields 


v 


= V, 


(C4) 


For  a 90°  magnetometer  (magnetometer  axis  perpendicular  to  rocket  axis),  the 
cos  ^ is  related  to  spin  and  coning  of  the  rocket  as  follows  (see  Figure  C2) 
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cos  <p  * cos  a cos  X » sin  « sin  X cos  uJat 

FOR  90*  MAGNETOMETER  X»90* 

AND  cos  ^ ■ sin  /cos  wjt 

Figure  C2.  Magnetometer  Geometry  for  Spinning  Rocket  in  a 
Magnetic  Field 


cos  </>  = sin  7 cos  u)gt  (C5) 

where 

y = instantaneous  angle  between  the  rocket  axis  and  the  earth's  magnetic  field 
ug  = spin  radian  frequency. 

Combining  Eq.  (C5)  and  (C4)  yields  the  instantaneously  voltage  output  of  the  mag- 
netometer in  terms  of  spin  and  coning 

ar  B 

v = Vfa  + sin  7 cos  ugt  . (C6) 

Since  sin  y (representative  of  coning)  is  a slowly  varying  function  in  time  com- 
pared to  cos  ugt  (the  spin  function),  Eq.  (C6)  is  an  amplitude  modulated  waveform 
(see  Figure  C3).  The  envelope  of  the  carrier  is  defined  by  two  equations  derived 
from  Eq.  (C6)  when  cos  wgt  = +1  and  -1  respectively 


0 


v = V. 
u b 


, AB  , B . 

+ -ir  + irsin  7 


v +*5. 

vb  k 


B 

"F'sin  7 


(C7) 

(C8) 
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Figure  C3.  Schematic  Showing  Typical  Magnetometer  TM  Volts 
Versus  Flight  Time,  = Upper  voltage  envelope,  vi  = lower 
voltage  envelope,  v=  instantaneous  magnetometer  output 


where 

vu  = upper  envelope  voltage 
Vj  = lower  envelope  voltage. 

The  average  value  of  the  envelope  voltage  at  any  instant  of  time  is 


V + V 

v =_y L-v  +*B 

ave  2 Vb+-TT 


Using  Eq.  (C9)  in  Eq.  (C6)  yields 


B 


v = v 


o . 


ave  + -k~sin  y cos  u t 


(C9) 


(CIO) 
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Likewise,  the  envelope  Eq.  (C7)  and  (C8)  become 


v 

u 


ave 


B 

+ “F'sin  7 


V1  =vave  ■~ITsin7  • {C11) 

Notice  that  vy  approaches  v^  only  when  y approaches  0 or  180°.  This  condition 
occurs  only  when  the  rocket  axis  is  aligned  with  the  magnetic  field.  Using  Eq. 

(C 11)  and  the  rocket  magnetometer  output  voltage  data  (vu  and  Vj),  the  rocket 
aspect  angle  (y)  relative  to  the  earth's  magnetic  field  can  be  found  from 


7 * 


k - 
2 B (vu 


(C12) 


In  the  northern  hemisphere,  the  angle  y is  usually  greater  than  90°  and  it  is  com- 
mon practice  to  describe  rocket  motion  in  terms  of  pitch  angle.  Pitch  angle  (P.A. ) 
is  defined  as  follows: 


P.A.  =180-7 


(C13) 


For  optical/ infrared  rocket  measurements,  it 
is  important  to  know  the  rocket  aspect  with  respect 
to  the  local  vertical  so  that  the  instrument  viewing 
angle  relative  to  the  vertical  can  be  determined. 
From  the  geometry  shown  in  Figure  C4  (spherical 
triangle  AEFA),  it  can  be  seen  that 

cos  7 = cos  6 sin  a cos  (D  - <P)  - sin  6 cos  a 

(C14) 

where 

7 * angle  between  rocket  axis  and  earth's 
magnetic  field 

a = angle  between  rocket  axis  and  the  local 
verticle 


6 * local  earth's  magnetic  field  dip  angle 
Figure  €4.  Spherical  (inclination) 

Geometry  of  Coning 
Rocket 
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D = local  earth's  magnetic  field  declination  angle 
0 = azimuth  of  rocket  axis  relative  to  true  north. 
Equation  (C12)  can  be  solved  for  sino  as  follows: 


sin  a = 


ab  ± c 


J72  2,2 

vb  - a + c 
b2  + c2 


(C15) 


where 

a * cos  7 

b = cos  6 cos  (D  - 0) 
c = sin  6 . 

Since  the  earth's  magnetic  field  parameters  6 and  D are  usually  known,  Eq.  (Cl5) 
can  be  solved  for  the  desired  rocket  aspect  angle  (a)  provided  the  rocket  heading 
(0)  is  known.  In  most  actual  cases,  the  rocket  heading  (0)  is  continually  changing 
with  time  as  the  rocket  precesses  (cones)  so  unless  other  data  (such  as  gyro  data) 
are  available,  it  is  difficult  to  prescribe  the  variations  in  rocket  heading.  A more 
practical  expression  for  a car  be  derived  in  terms  of  the  terms  of  the  rocket  axis  of 
precession  since  the  heading  of  this  axis  is  not  as  susceptible  to  change. 

To  derive  this  relationship,  the  angle  a is  defined  in  terms  of  the  angle  (6^) 
between  the  magnetic  field  and  the  rocket  axis  of  precession  and  the  coning  half 
angle  (p).  If  it  is  assumed  that  the  rocket  coning  is  sinusoidal,  then  from  the 
geometry  of  Figure  C4  (spherical  triangle  AEDA) 

cos  7 = cos  <?c  cos  P + sin  0 Q sin  p cos  (u>ct  - U ) (C16) 

where 

7 = angle  between  rocket  axis  and  earth's  magnetic  field 

0^  = axis  of  precession  (coning)  relative  to  the  earth's  magnetic  field 

p = half  cone  angle 

«c  = coning  radian  irequency 

t * time 

fi  = phase  angle. 

The  maximum  and  minimum  values  of  7 can  be  found  from  Eq.  (C 16)  when 
cos  (wfit  - fi)  = -1  and  +1  respectively,  thus 
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(C 17 ) 


y . = 0 - fj 

'min  c ' 

Both  (3  and  $c  can  be  found  from  knowledge  of  y and  Eq.  (C17)  as  follows: 


7 - y 

- max  'min 


P - - 


7 + v 

n 'max  'min 
9c g * 

From  geometry  in  Figure  C4  (see  spherical  triangle 

cos  0 = cos  6 sin  a cos  (D  - 0 ) - sin  6 cos  a 
c c c c 

which  can  be  solved  for  sin  a 


sin* 

c 


„v 


n,2  -l2+n2 


m 


2 +n2 


ADFA), 


(C 18) 


(C 19) 


(C20) 


where 

1 = cos  0 

c 

m = cos  6 cos  (D  - ^c) 
n = sin  6 

0c  = azimuth  of  axis  of  precessions  relative  to  true  north 
«c  = angle  of  axis  of  precession  relative  to  local  vertical. 

All  other  symbols  have  been  defined  previously. 

The  phase  angle  (Cl)  in  Eq.  (C16)  and  shown  in  Figure  C4  can  be  described  in 
terms  of  the  heading  (^c)  of  the  precession  axis  as  follows.  Referring  to  Figure 
C4,  the  spherical  right  triangle  ADFA  yields 


sin  (90  +6)  i; 
sin  (180  - Cl) 


sin  0 
c_ 

sin 
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which  reduces  to 


sin  H = 


sin  A0c  cos  6 
sin  6 


(C22) 


where 


A0c  = D - ^ 


All  other  symbols  have  been  defined  previously.  The  phase  angle  can  also  be 
derived  in  terms  of  a , 6c  and  where 


cos  f2  = ■ 


cos  6 cos  A0  - cos  6 sin  a 


sin  0 cos  a 
c c 


(C23) 


The  aspect  anjle  (p)  of  the  rocket  axis  relative  to  the  local  vertical  is  (see 
Figure  C4)  spherical  triangle  EDFE  where 


cos  a = cos  a cos  8 - sin  a sin  P cos  w t . 

c c c 


(C24) 


In  many  actual  cases,  the  half  cone  angle  (£)  is  small  such  that  the  following 
approximations  are  possible: 

cos  /3  ~ 1 


sin  p ~ P 


(C25) 


For  this  case,  Eq.  (C24)  and  (CIS)  become 


cos  a ~ cos  a c - p sin  »c  cos  i^,t 


(C2S) 


or 


a ~ a +8  cos  w t 
c c 


cos  y ~ cos  0 +8  sin  9 cos  (w  t - Q) 
c c c 


(C27) 

(C28) 


or 


y ~ 9 - p cos  (w  t - il) 
c c 


(C29) 
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If  the  additional  approximation  can  be  made  that  the  phase  angle  (57)  is  small  and 
can  be  neglected  (see  Eq.  C23),  then  0 cos  uct  can  be  eliminated  from  Eq.  (C27) 
and  (C29)  to  yield 

a ~ c + - y . (C30) 

The  error  (e)  encountered  by  using  the  approximation  (Eq.  (C30))  can  be 
derived  by  adding  Eq.  (C27)  and  (C29)  and  making  use  of  the  identity 

cos  uct  - cos  (uct  - SI)  * -2  sin  (57/2)  sin  (uct  - 52/2)  . 

The  instantaneous  error  (e)  is 

c = -2/3  sin  (0/2)  sin  (u  t - 0/2)  (031) 

and  thus  the  maximum  error  (e„.„)  is 

m six 

emax=±2^sinf  • (C32) 

For  many  cases  of  well  behaved  free  ballistic  rocket  flight,  the  angles  defining 
the  rocket  axis  of  precession,  ac  and  9C>  are  constant  for  at  least  some  portions  of 
flight.  For  these  cases,  a convenient  expression  can  be  found  for  rocket  aspect  a 
using  Eq.  (C30)  and  (C13): 

a ~ P.A.  + Constant  (C33) 

where 

Constant  = or  + 0 - 180°  . 

c c 

Equation  (C32)  is  particularly  useful  since  the  maximum  angular  error  can  be 
estimated  by  Eq.  (C32).  The  procedure  for  evaluating  Eq,  (C33)  follows: 

(1)  Compute  the  angle  y and  P.A.  using  the  magnetometer  data  and  Eqs.  (Cl2) 
and  (C13). 

(2)  Compute  9 and  0 from  the  time  history  of  y found  in  step  (1)  using 
Eq.  (C18). 

(3)  Compute  »c  from  9 ^ and  a known  or  estimated  value  of  0c  (heading  of  the 
axis  of  precession  using  Eq.  (C20). 

(4)  Compute  the  estimated  aspect  angle  a from  a and  0 and  the  P.A.  data 
found  in  (1)  above  using  the  approximation  Eq.  (C33).  The  maximum  angular  error 
in  or  can  be  found  by  evaluating  Eq.  (C22)  and  (C32). 
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It  must  be  pointed  out  that  there  is  an  additional  error  in  aspect  angle  a not 
accounted  for  by  Eq.  (C31)  and  this  arises  from  the  assumption  that  the  cone  angle 
(3  is  small.  This  can  be  illustrated  by  an  example  where  the  heading  of  the  rocket 
axis  of  precession  coincides  the  heading  of  earth  magnetic  field,  that  is  0C  = D. 

For  this  case 

a + 6 = 90°  + 6 
c c 

and  Eq.  (C33)  reduces  to 
a = P.A.  + 6 - 90°  , 

where  6 is  the  magnetic  field  dip  angle.  Also,  for  this  case  Eq.  (C22)  indicates 
that  Cl  = 0 and  thus  Eq.  (C32)  would  yield  a maximum  error  of  0°.  This  result  is 
obviously  incorrect  but  for  cases  where  (3  is  small  (£6°),  the  discrepancy  can  be 
shown  to  be  negligible  for  most  applications.  For  instance,  assume  a measured 
P.A.  to  be  50°  at  WSMR  (6  = 60,  D * 10)  and  (3  = 6°.  The  approximation  fora  yields 
a = 20~.  On  the  other  hand,  Eq.  (C15)  can  be  solved  fora  without  approximation. 
For  the  worst  case  the  instantaneous  rocket  axis  heading  (0),  as  distinguished  from 
the  axis  of  precession  heading  (0c>,  differs  from  the  magnetic  field  heading  (D)  by 
the  half  cone  angle  (fi );  that  is 

D ±/3  = 0 . 

When  Eq.  (C15)  is  evaluated  under  tht.,e  conditions,  the  value  of  a is  20.  07°  - a 
discrepancy  of  only  0.  07°. 
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Appendix  D 

Derivative  Using  Name  fieri  Filtering* 


A numerical  filter  generally  consists  of  a set  of  weights  w{n)  applied  to  a set 
of  input  data  (see  I below)  points  XU)  to  form  a "filtered"  set  of  data  Y (I). 

For  this  problem,  an  imperical  number  N was  set  and  a low  pass -derivative 
filter  was  calculated  and  applied  to  the  data  as  follows: 

N 

w(n)X(n  + f)  . I 

n-N 


The  values  of  w(n)  were  calculated  as  follows:  An  array  of  8192  values  of  frequency 
response  versus  frequency  W(f)  was  set  up  on  the  computer.  The  array  consisted 
of  8192  complex  values  at  equally  spaced  data  values  from  -4095  Af  to  4096  Af 
where 


* ' 8192  At 

and  At  is  the  tune  spacing  of  the  input  data.  The  value  at  each  frequency  (f)  was 
-(2*f)i  C(f),  where  -(2*f)i  represents  the  time  (t)  derivative  of  e 1 divided  by 

e(~2*f)i  t and  C(f)  is  the  response  curve  used.  For  this  procedure: 

*Thi8  analysis  has  been  accomolished  by  Kiel  Grossbard  of  Boston  College  for 
AFCRL  Data  Computer  Center. 
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C(f)  = 1 


for 


If!  sf, 


m.n 


= l/2  cos  — 


2*<lfl  - If  J> 


min’ 
max  ^min^ 


+ 1 for 


f . < |f|  < f 
mm  1 1 max 


= 1 


for  |f!  2:  f 


max 


where  f and  f are  imperical  input  frequencies, 
min  max  ^ e 

The  time  series  representation  of  W(f),  (U(t))  was  then  derived  using  a numer- 
ical Fast  Fourier  Transform.  The  results  for  U(-N  At)  to  U(+N  At)  formed  the 
w(n)  between  -N  and  +N  used  in  the  computer  program  to  calculate  diQ/dt. 

The  resultant  filtering  only  approximates  the  frequency  response  function  C(f) 
for  two  reasons.  The  first  reason  is  the  truncation  of  U(t)  such  that  W(n)  = 0 for 
I n | > N.  This  "error"  minimized  by  increasing  N until  the  final  results  of  the  pro- 
gram did  not  change  by  1 percent  with  increasing  N.  The  second  reason  is  that  the 
function  is  only  set  at  the  values  n Af  for  n = -4095  to  4096  and  not  at  frequency 
values  falling  between  these  points.  This  problem  was  minimized  by  using  many 
values  (8192)  in  W(f)  and  by  using  a "smooth"  C(f).  In  particular,  the  behavior  of 
C(f)  between  fmin  and  f was  set  up  to  form  a "smooth"  change  from  a derivative 
response  to  a zero  response.  This  error  was  not  checked  but  experience  indicates 
that  this  error  should  be  less  than  3 percent  for  the  parameters  used.  In  any  case, 
the  errors  due  to  the  inaccuracy  in  fitting  C(f)  are  probably  less  than  the  error  in 
assuming  C(f)  specifies  the  necessary  frequency  response  curve  for  reducing  X(t). 

Using  this  analysis,  we  find  dRQ/dt.  These  results  were  transformed  to 
dRQ/dh  as  follows: 

Let  Z(/)  be  the  series  of  values  of  dR  / dt  with  corresponding  time  array  tt(n) 


and  height  array  h(i  ) and  V jj  array.  V^U ). 


N (h(i))  = - ^-r) 
2 eff  ' 


Z(f ) 


tt(f  + 1)  - tt(f  - 1) 
hU  + l)  - h(/  - 1) 


where 


ttfi  + 1)  - tt(f  - 1) 
h(t  + 1)  - h(jf  - 1) 
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is  used  to  approximate  ^ and  the  chain  rule 

dR_  _ dR  dt 
dh  ~ dt  dh 

has  been  used.  This  approximation  for  dt/dh  is  another  small  source  of  error  in 
the  values  of  Nz(h(f )).  This  error  should  be  minimized  by  the  slow  change  of  h{i) 
and  tt(f ) in  the  input  data. 

In  order  to  make  an  estimation  of  the  error  in  M z (h ) > » the  following  proce- 
dure was  used: 

To  find  the  error  for  N (h(f))  consider 
z 

N (h(f  - D)  + N (h(/ ))  + N (h(l  + 1) 

S 5 " NzM<h<*» 


i + 1 


ctn  (h(f ))  = 


E (Nz(h(j))-NzM(h(f))2 

jg':*  . 


Then  it  is  assumed  the  correct  answer  for  Nzc(h(/))  (NzM(h(i)))  lies  in  the  region 


zM' 


NzM(h(i))  ' aNzM(h(i))s  Nzc<h</,)  ~ NzM{h(i))  +aNzM<h(i>) 


with  a probability  -M3.  68. 

This  procedure  should  give  an  overestimate  of  the  error  if  the  results  of 

N (hU))  are  unbiased  since  any  change  between  N (h(f  - 1)),  N (h(f ))  and 
z z c z c 

N (h(f  + 1))  tends  to  increase  calculated  value  ofow  (h(i )).  At  the  turning  points 

of  the  N (h(f )),  however,  this  method  tends  to  underestimate  the  strength  of  the 
zc 

peak. 
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